arXiv: 1501.0509lv4 [physics.soc-ph] 1 Sep 2015 


Mapping the (/-voter model: From a single chain to complex 

networks 


Arkadiusz J§drzejewski a , Katarzyna Sznajd-Weron a , Janusz Szwabinski b 

a Department of Theoretical Physics, Wroclaw University of Technology, Poland 
b Institute of Theoretical Physics, University of Wroclaw, Poland 


Abstract 

We propose and compare six different ways of mapping the modified q-voter model to 
complex networks. Considering square lattices, Barabasi-Albert, Watts-Strogatz and 
real Twitter networks, we ask the question if always a particular choice of the group of 
influence of a fixed size q leads to different behavior at the macroscopic level. Using 
Monte Carlo simulations we show that the answer depends on the relative average path 
length of the network and for real-life topologies the differences between the considered 
mappings may be negligible. 

Keywords: Opinion formation, Opinion dynamics, Q-voter model, Agent-based 
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1. Introduction 

Ordering dynamics is not only a classical subject of non-equilibrium statistical physics 
0,1, but also one of the most studied issues in the field of sociophysics lii- It 
often represents opinion dynamics under the most common type of the social influence, 
known as conformity. Among many others, models with binary opinions are of particular 
interest. 

One of the most general models of binary opinion dynamics was introduced by Castel¬ 
lano et al. @ under the name the q-voter model as a simple generalization of the original 
voter model jJ. In the proposed model, q randomly picked neighbors (with possible 
repetitions) influence a voter to change its opinion. If all q neighbors agree, the voter 
takes their opinion; if they do not have an unanimous opinion, the voter can still flip 
with probability e. It has been argued that for q = 2 and e = 0 the q-voter model co¬ 
incides with the modified Sznajd model, in which unanimous pair (in one dimension) of 
the neighboring sites SiSi+i influences one of two randomly chosen neighbors i.e. S)_i 
or S i+2 0. 

Following this reasoning in Ref. Q we have introduced a modified one-dimensional 
version of the q -voter model, as a natural extension of the original voter and the Sznajd 
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model: a panel of q neighboring spins Si, Si+i, Si+ q _i is picked at random. If all 
q neighbors are in the same state, they influence one of two neighboring spins Si-\ or 
Si+q- If not all spins in the q-panel are equal then nothing changes in the system. This 
modification has been later considered in Refs. 


For e = 0 both formulations of the q -voter model seem to be almost identical with the 
exception of the repetitions possible in the original version [&!]. However, there is another 
difference between formulations [b] and :8|, namely the first belongs to the class of so 
called inflow and the second to outflow dynamics BOH- There was a controversy 
related to the subject if the inflow and outflow dynamics are equivalent [8. 13. 15. lil 171. 
Recently, it has been shown that they are equivalent for q = 2 0 , at least in respect to 
the exit probability. However for larger values of q, even in one dimension the situation 
is not clear ITcj [ill. |T3. ITij . Moreover, differences between dynamics in respect to the 
phase transitions induced by the stochastic noise has not been investigated up till now. 

At first glance, it seems to be trivial that choosing a different set of interaction 
partners will lead to different results on the macroscopic scale. However, as described 
above, the problem occurred to be not as simple as it seems and gained the attention in 
the literature. Therefore, one of the aims of this paper is to contribute to the outflow- 
inflow discussion. 

The second, probably more significant, aim is related to applications of the q-voter 
model in social sciences, because as noted by Macy and Wilier there was a little effort to 
provide analysis of how results differ depending on the model designs 0. Moreover, in 
respect to social applications one could ask the question - how to construct the group 
of influence to create easier an order (consensus) in the system? Another question one 
could ask - is the problem relevant for any type of a network or maybe for some networks 
different types of the group of influence will lead to the same results at the macroscopic 
scale? Only in the case of a complete graph the definition of the q-voter model is straight¬ 
forward - since on this topology all spins are neighbors 2C)|, 211, all proposed versions of 
the q -voter model are equivalent. 

In the context of opinion dynamics it would be however desirable to consider the 
models on top of more complex networks, as they are better representations of contact 
patterns observed in the social systems [22|, [23| • There are already several attempts to 

However, as shown in 
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generalize the g-voter model to complex networks 

Ref- [13, even in the simple case of transferring the model from ID chain to a 2D square 
lattice there is no unique rule of choosing the group of influence. Thus, the main goal 
of this paper is to check, how different ways of picking up the group may 
impact the macroscopic behavior of the model. Specifically, we will focus here 
on the phase transitions induced by the stochastic noise that represents one type of the 
social response, known as independence dJSij. 


2. Model 

Within the modified q-voter model we consider a set of N agents called spinsons. 
This name, being a combination of the words “spin” and “person”, is used to emphasize 
that the Ising spins in our model represent persons characterized by only one binary trait 
(a detailed explanation of this notion may be found in Ref. EH).' Each i-th spinson has 
an opinion on some issue that at any given time can take one of two values S t = ±1, i = 
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1,2,..., JV (“up” and “down”). The opinion of a spinson may be changed under the 
influence of its neighbors according to two different types of the social response [28| : 


• Independence is a particular type of non-conformity. It should be understood as 
unwillingness to yield to the group pressure. Independence introduces indetermi¬ 
nation in the system through an autonomous behavior of the spinsons [29| . 


• Conformity is the act of matching spinson’s opinion to a group norm. The nature of 
this interaction is motivated by the psychological observations of the social impact 
dating back to Asch [30|: if a group of spinson’s neighbors unanimously shares an 
opinion, the spinson will also accept it. 


Other types of the social response are possible as well (see Ref. [2l| for an overview), but 
the above two are of particular interest for studying opinion dynamics. 

We study the model by means of Monte Carlo simulations with a random sequential 
updating scheme. Each Monte Carlo step consists of N elementary events, each of 
which may be divided into the following steps: (1) pick a spinson at random, (2) decide 
with probability p , if the spinson will act as independent, (3) if independent, change its 
opinion with probability 1/2, (4) if not independent, let the spinson take the opinion of 
its randomly chosen group of influence, provided the group is unanimous. More details 
on the dynamic rules of the model may be found in Ref. [2(J. It is worth to stress here the 
difference between the modified g-voter model with independence p > 0 and the original 
q -voter model with e > 0 jfjj] . In fact one could introduce a generalized model with both 
parameters p > 0 and e > 0, in which each elementary time step is described by the 
following algorithm: 


1. Choose at random one spinson Si located at site i. 

2. Decide with probability p, if the spinson will act independently. 

3. In case of independence, a spinson flips to the opposite state with probability 1/2. 

4. In other case (conformity), choose q neighbors of site i (a so called g-panel): 

(a) If all the q neighbors are in the same state, i.e. g-panel is unanimous, the 
spinson takes the state of the q neighbors. 

(b) Otherwise, i.e. if g-panel is not unanimous, spinson flips with probability e. 

Clearly, the original g-voter model is a special case of the above algorithm with p = 0 
and the model considered here is a special case with e = 0. Note, that in contrast to p, 
parameter e does not describe the independence. For p = 0, if only the unanimous g-panel 
exists, the spinson will take its state, which means that it never acts independently. In 
consequence, the state with all spinsons in the same state is the absorbing steady state 
for the original g-voter model, whereas it is not for the model considered here unless 
p = 0. Therefore, the original g-voter model with p = 0 is not suitable to model e.g. 
diffusion of innovation, for which the initial state with all spinsons down (unadopted) is 
a typical one 0 . 

Introduction of a group pressure as one of the rules governing the dynamics assumes 
some form of interactions between the spinsons. Those interactions are best illustrated 
as connections between nodes of a graph the spinsons are living on. In its original for¬ 
mulation, the g-voter model can be easily investigated on an arbitrary lattice and the 
definition is clear, since q individuals influencing the voter are chosen with repetitions 
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from the nearest neighborhood of the voter. Therefore, even in one-dimension the param¬ 
eter q can have an arbitrary value [(|. Although such a definition of the model leads to 
interesting results from the physical point of view, it seems to be sociologically unreliable. 
In the modified version, repetitions are forbidden, which is probably more sociologically 
justified. Moreover, q influencing agents may form panels of different kinds, including 
structures proposed in fl2l . [26}, which is also different from the original formulation of 
the model. This in turn allows to investigate the role of the group structure. However, 
such a modified definition causes ambiguity in mapping the model on an arbitrary graph. 

We use here both the Watts-Strogatz [31]] and the Barabasi-Albert networks [32| as 
the underlying topology of spinson-spinson interactions, since they nicely recover the 
small world property of many real social systems [33}. We set q = 4 for two reasons: 
(1) to reflect the empirically observed fact that a group of four individuals sharing the 
same opinion has a high chance to ’convince’ the fifth, even if no rational arguments 
are available [3(3, [34j and (2) to compare our results with those obtained on the square 
lattice. 

Chosen from the plethora of possibilities, in Fig. [T] six different groups of influence 
on a complex network are schematically shown. We would like to stress here that the 
choice of precisely such groups is not accidental and is dictated mainly by earlier papers 
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Line - after picking up a random target spinson (marked with a double red circle 
in the figure), we randomly choose one of its neighbors, then one of the neighbors 
of the neighbor and finally a neighbor of the latter one. All members of the g-panel 
are indicated with a blue circle in the figure. This is the natural generalization of 
the ID g-voter model and was used e.g. in Ref. 
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• Block - the group consists of a random neighbor of the target spinson, and three 
neighbors of the neighbor. This method resembles to some extent the 2 x 2-block 
used on square lattices in 2D and was used for instance in Ref. '26j. 


• NN - four randomly chosen nearest neighbors of the target spinson are in the 
group. This method was used in the original g-voter model 


• NN3 - this is a slight modification of the NN method leading to an extended range 
of the influence: the group is composed of three randomly chosen nearest neighbors 
of the target spinson, and a neighbor of one of those nearest neighbors. We have 
introduced this method just to investigate the impact of the range of an influence 
group on system’s ability to stay in an ordered state. 


• RandBlock - a spinson and its three neighbors build the the group of influence 
as in the Block method. However, the block may be located anywhere on the 
network. This method has been chosen as a reference for the mean-field type 
approach represented by the next method and similarly as NN 3 it is aimed to 
investigate the role of the range of interaction. 


• Rand - the group consists of four randomly chosen spinsons, not necessarily con¬ 
nected with the target spinson. This corresponds to the mean-field app roach, for 
which analytical results on the phase transitions are already known [2(1 ]. 
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Figure 1: (Color online) Different groups of influence on a complex network. A random spinson is marked 
with a double red circle, the spinsons in its the group of influence - with a single blue one. Interactions 
within the group and between the target spinson and the group are represented by thick lines. See text 
for explanations. 
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Figure 2: Magnetization m as a function of the independence factor p for different groups of influence 
on the square lattice N = 100 X 100. As seen, there is no phase transition for NN model on the square 
lattice and the critical value of p increases with the interaction range, as expected. 


Note that in case of RandBlock and Rand methods we actually abstract away from 
the underlying network topology of the model. We expect both methods to be equivalent 
to the complete graph case if the minimum degree of a node in the network is bigger 
than or equal to q = 4 (otherwise the RandBlock may differ slightly from the complete 
graph, because there will be not always enough spinsons in the neighborhood to build the 
influence group). Although we use Rand mostly as a benchmark for our simulations, the 
RandBlock is much more interesting because it corresponds to a situation often encoun¬ 
tered in many organizations, in which an informal and unknown network of interactions 
is over imposed on the given formal communication structure. 


3. Results 


The main goal of this paper is to answer the questions if and how details at the 
microscopic level manifest at the macroscopic scale. Among other macroscopic phenom¬ 
ena, phase transitions are certainly the most interesting ones. For the models of opinion 
dynamics, the most natural order parameter is an average opinion to, defined as magne¬ 
tization i.e. to = jf y~) Si . It has been shown that in the case of theo-voter model the 
phase transition may be induced by the independence factor p |29l 120 ]. Below the crit¬ 
ical value of independence, p < p c , the order parameter to ^ 0. For high independence, 
p > p c , there is a status-quo, i.e. to = 0. Such results were obtained on the complete 
graph topology which corresponds to the mean field approach, as well as on the square 
lattice but only for the one particular choice of the q-panel equivalent to the Sznajd 
model [29]. In [ 29 } a 2 x 2 box of four neighboring spinsons were chosen randomly and 
influenced one of the 8 neighboring sites of the box. Here we test six different methods 
described in the previous section (see Fig. [2]). 

It is seen that the critical value of the independence factor strongly depends on 
procedure of choosing an influence group. The phase transition is observed for all meth¬ 
ods except of NN, which corresponds to so called inflow dynamics 14J. As expected, 
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Figure 3: Magnetization m as a function of the independence factor p for six different groups of influence 
on the Barabasi-Albert network of size N = 10 4 and parameters Mo = 8 and M — 8. Four models (Block, 
Line, RandBlock and Rand) collapse into a single curve and only two (NN and NN3 models) can be 
distinguished from others. 


RandBlock and Block methods overlap and agree with MFA result found in [20[, i.e. 
p c (q) = (<Z — l)/(g — 1 + 2 9_1 ), which for q = 4 gives p c = 3/11. Methods for which the 
range of interaction is shorter tend to show lower critical value of p. This result is very 
intuitive, since the infinite range of interactions usually corresponds to MFA and gives 
the largest critical value. 

Results on the Barabasi-Albert (BA) network (see Fig. [3J are less intuitive. It 
occurs that for this topology differences between methods are almost negligible. The 
phase transition is observed for all six models and the critical value of p changes only 
slightly with method. Four models (Block, Line, RandBlock and Rand) collapse into a 
single curve and only two (NN and NN3 models) can be distinguished from others. The 
natural question arises here - why differences between models are clearly visible on the 
square lattice and are almost negligible in the case of BA? 

It should be recalled that the average path length l for the square lattice increases 
with the system size TV as l ~ TV 1 / 2 , whereas in the case of BA as l ~ InN/lnlnN 
{23 |. It means that for the same system size the average path length is dramatically 
shorter on BA than on the square lattice. In result the range of interactions on BA is 
effectively much larger. To check the role of the average path length we have simulated 
all 6 methods on the Watts-Strogatz network. This topology is particularly convenient 
because for the fixed network size TV it is possible to decrease the average path length l 
by increasing the rewiring probability /?. 

Results for several values of /? are presented in Fig. [4] As /? increases the critical 
point p c shifts towards higher values and simultaneously differences between 4 methods 
vanish up to a threshold value f3 = 0.5. Results for larger values of (3 (i.e. (3 > 0.5; not 
shown in Fig. II are identical with those obtained for /3 = 0.5. To check how results 
scale with the system size we simulated models on networks of sizes from TV = 500 to 
TV = 10 5 (see Fig. [5J. Surprisingly, results are virtually independent on the system size. 
Analogous results has been obtained for the exit probability in one dimensional system 
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Figure 4: Magnetization m as a function of the independence factor p for different groups of influence 
on the Watts-Strogatz network of size N = 10 4 with the average degree K = 8 and rewiring parameter 
(3. The critical value of independence p c increases with the range of interactions, as expected. However, 
with increasing (3 differences between models vanish and for large [3 only two models (NN and NN3) can 
be distinguished from others. 


with inflow HI and outflow dynamics fid . 

The fact that results do not depend on the system size undermines our predictions 
that the average path length l itself determines if all mapping methods overlap or not, 
because the l increases with the system size [35l ]. However, one should probably not look 


at the path length itself but at the relative path length, which is defined as the average 
path length of a given network divided by the average path length of a random network 
of the same size and average degree. Normalizing networks’ characteristics by those of 
the correspond ing random graphs is a procedure usually used to compare networks of 
different sizes [36, 37, 


Thus, we will use the relative path length to describe the 


networks under consideration and to compare them. For example for the Watts-Strogatz 
network with k = 8 and /3 = 0.05 the relative path length is equal l re i = 1.65422863636 
for N = 500 and l re i = 1.67766421808 for N = 1000, i.e. almost size independent. For 
Barabasi-Albert of size N = 500 it is much shorter i.e. l re i = 0.974 and almost does not 
change with the system size. Interestingly, if one consider relative path length l re i for 
Watts-Strogatz of size N = 500 and k = 8 with different /? it occurs that relative path 
length approaches 1 for /3 = 0.5 and precisely: l re i{P = 0.05) = 1.65422863636, l re i{P = 
0.1) = 1.39668726543, l re 0 = 0.2) = 1.2034321044, l re 0 = 0.5) = 1.02497507194. This 
somehow explains why results for BA and WS with /3 = 0.5 are almost identical. 

To check if our predictions about the relative path length are correct, we decided 
to investigate the problem on several networks of the same size, including real Twitter 
networks. We took Twitter data from the Stanford Large Network Dataset Collection 
available at https://snap.stanford.edu/data/egonets-Twitter.html [39|, because 
it includes about 1000 different networks with a broad spectrum of diverse characteristics. 
It was relatively easy to find in the dataset networks of the same size, but with different 
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Figure 5: Magnetization m as a function of the independence factor p on the Watts-Strogatz network 
with the average degree K = 8 and rewiring parameter /3 = 0.05 for several system sizes. Results are 
not influenced significantly by the system size and the critical value of independence p c increases with 
the range of interactions, as expected. 


average path lengths and/or clustering coefficients. Thus the dataset was well suited for 
testing our hypothesis about the path lengths. Magnetization m as a function of the 
independence factor p on six different networks of size N = 233 is presented in Fig. [6] 
In the top row results on three real Twitter networks [39J are presented. Networks in the 
left and middle top panels have almost identical path length l but different clustering 
coefficient C. On the other hand, middle and right networks have almost identical 
clustering coefficient but slightly different path length. It seems that results for all 
methods overlaps the best on the right network, which has the shortest path length. 
Simultaneously, it seems that results on left and middle networks are the most similar 
to each other, i.e. path length l is more significant than clustering coefficient C in 
determining if all mapping methods will give the same result or not. However, because 
the differences between properties that we take into account (i.e. I and C) do not vary 
much from network to network, for all three Twitter networks almost all methods collapse 
into a single curve and only NN and NN 3 slightly deviate from others. In the bottom 
row of Fig. [0] results on three artificial networks are presented. We took two Watts- 
Strogatz networks (left and middle plots) and a Barabasi-Albert one. In case of the WS 
network with the clustering coefficient similar to Twitter networks and much longer path 
length l, each mapping gives a completely different result (bottom left). On the other 
hand, if the path lengths are similar to those of real networks, for both WS and BA 
models we observe already known behavior - most methods collapse into one curve and 
only NN and NN 3 differ slightly from the others. Thus, this feature does not depend on 
the network topology and the clustering coefficient, which again confirms that the path 
length is the significant property for the investigated problem. 

From obtained results we also conclude that nor the average degree neither the degree 
distribution are significant for the investigated problem. We realize that there are more 
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Figure 6: Magnetization m as a function of the independence factor p on six different networks 
of size N = 233. In the top row results on three real Twitter ego networks f39tl are presented 
whereas in the bottom row results on three artificial networks are shown. The Twitter network IDs 
in the titles of the plots correspond to file names (‘ego’ node IDs) in the Twitter dataset taken from 
https://snap.stanford.edu/data/egonets-Twitter.html [39tl . In the bottom row results on three 
model networks are shown: a Watts-Strogatz network with the average degree K = 4 and the rewiring 
probability (3 = 0.035 (bottom left), a WS network with K = 18 and (3 = 1 (bottom middle) and a 
Barabasi-Albert one with the number of new edges M = Mq = 10 (bottom right). 
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properties of the real networks that could be taken into account, but our prediction about 
the importance of the relative path length seems to be reasonable also from the statistical 
physics point of view. It should be noted that for all networks considered in this paper 
two methods collapse into a single curve - RandBlock and Rand and overlap mean-field 
results presented in [20| ■ This can be easily understood in the case of Rand method - 
regardless of the topology, neighbors are randomly chosen. In the case of RandBlock, 
obtained results are less obvious. However, this method introduces interactions with the 
infinite range. We know from statistical physics that the mean field approach should 
give exact results in the case of infinite interactions. Following this reasoning we can also 
understand why, with decreasing path length, results for all mapping methods approaches 
the mean field results - the relative interaction range increases. Another phenomena that 
can be understood on this basis is the fact that TV TV method differs the most from the 
Rand and Line is the most similar. TVTV method has relatively the shortest range of 
interactions and Line much larger. 

The differences between methods may be also explained, at least qualitatively, in 
terms of probabilities of finding non-unanimous influence groups. For the sake of sim¬ 
plicity let us assume that our network has the topology of a Bethe lattice [4(j] with the 
coordination number fc. For convenience, we took a modified definition of the Bethe 
lattice with the central node having only k — 1 neighbors (Fig. [7]). Thus, the central node 
has k — 1 neighbors in its closest neighborhood, ( k — l) 2 agents at the second level of its 
ego graph, and in general (k — l) d nodes at distance d. Now, let us consider our model at 
an early stage of a simulation. Let us assume that there are only two spinsons including 
the central one in the “down” state due to independence and that the central spinson has 
been chosen again in a basic Monte Carlo event (it will be referred as the target spinson 
in the following). However this time it is not independent, i.e. it is exposed to the group 
pressure. Since most of the spinsons are in the “up” state, the system has a natural 
tendency to reduce disorder due to conformity. Nevertheless, we can ask the question 
whether there are significant differences between the methods in maintaining disorder 
in the system. In other words we can check if the methods differ in the probabilities of 
finding a non-unanimous group of influence in this situation. 

To this end, we can consider configurations with the other “down” spinson residing 
at different levels of the target’s ego graph. We start with the “down” spinson being in 
the nearest neighborhood of the target node. In this case we expect the NN method to 
give the highest probability to build a group the “down” spinson belongs to. The reason 
is simple: this is the only method which operates exclusively in the nearest neighborhood 
of the target spinson. Thus we draw 4 agents out of fc — 1 to form the group and it is 
very likely that the “down” spinson will belong to the group (see top left plot of Fig. [7] 
for a schematic representation). The probability of finding a non-unanimous group is 
slightly smaller in TVTV3 case, because here only 3 drawings out of fc — 1 from the closest 
neighborhood are allowed. The Line and Block methods require only one drawing from 
the first level (top right plot of Fig. Q. Hence it is less likely to hit the “down” spinson. 
Finally, RandBlock and Rand algorithms yield the smallest probabilities, because they 
operate on the whole network rather than in the close neighborhood of the targ et node. 

Since the problem at hand is nothing but a variation of an urn problem [411 ]. we can 
actually calculate for each method the probability of finding a non-unanimous group of 
influence. To focus our attention we set fc = 9 and the system size TV = 4680 meaning 
that the ego graph of the target spinson consists of 4 levels. In the case of the NN 
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method the number of all possible g-panels is just the number of 4-combinations selected 
from the nearest neighborhood of the target, 



(1) 


The not adopted agent at the first level has to belong to each non-unanimous group. The 
other three members are selected from seven “up” spinsons residing at that level. Thus, 
the number of non-unanimous groups in the NN method is equal to 



( 2 ) 


This yields the following probability of finding a non-unanimous group: 


p(i) _ 
^NN ~ 


\NU\ 


= 0.5. 


( 3 ) 


The superscript (1) in the last expression indicates the level the other “down” spinson 
belongs to. In the NN3 method we select three agents from the first level, then pick one 
out of them and draw one spinson from its neighbors at the next level. This gives 





= 1344 


( 4 ) 


possible influence groups. Once the “down” spinson is chosen, we select two others from 
the first level, then pick one of the members of the group and add one of its neighbors 
to the group. The number of all possibilities is in this case given by 



Hence, 


p(i) 

*NN3 


\m 

|Q| 


0.375. 


Similar analysis leads to the following results for the other methods: 


( 5 ) 

( 6 ) 


Pmock = 0-125, L = 0.125, 

PrL = 0-85 * IQ" 3 , P^ ndB iock = 0-21 * 10- 3 . (7) 


We see that indeed the NN method gives the highest probability of leaving the target 
spinson untouched in this case. 

If the other “down” spinson resides at the second level of target’s ego graph, the 
Block method should give the highest probability of maintaining disorder in the system, 
because it consists of drawings mostly from that level (bottom row of Fig. Q. Again, we 
can calculate the corresponding probabilities to get: 

PbL = 0-046, P£ e = 0-015, P^i 3 = 0.015, 

4 2 <L = 0.85*10- 3 , Pg>ndBiock= 0-29 *1CT 3 , p™= 0. (8) 
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The Block method gives indeed the highest probability, followed by Line and NN 3. 
The probability for NN is zero, because the method operates only at the first level. It is 
interesting to note that the probabilities for the second level are much lower than those 
for the first one. With similar reasoning we can show that in general the farther the 
distance between the two “down” spinsons, the smaller the chance to maintain disorder, 
i.e. to let the target spinson unchanged. As a consequence, only the state of the two 
closest levels is actually significant for the evolution of target’s opinion. For this reason 
the methods NN and NN3, followed by Block and Line destroy the order in the systems 
the fastest, i.e. for relative small values of the independence parameter p. 

Note that the above conclusion is in accordance with our simulation results shown 
in Figs. [3j [I] and [6] Thus, although a Bethe lattice resembles ego graphs of model and 
real networks on average only, the reasoning remains the same - if the “down” agents are 
sparse, the NN and NN 3 methods yield the highest probability to maintain disorder in 
the system. 

From our simulations it follows that in networks with short relative average path 
lengths the differences between all but NN and NN 3 are negligible. A short path length 
means usually that the ego networks of all agents are rather “flat”, i.e. most agents reside 
at very few levels. In this case already the second and next levels of an ego graph are 
highly populated, leading to negligible probabilities of finding non-unanimous groups at 
the beginning of a simulation for methods operating mainly beyond the first level. Hence, 
on such networks the results for four methods (Block, Line , RandBlock and Rand) are 
essentially indistinguishable. The NN and NN 3 methods deviate slightly from others 
giving a bit smaller critical values of the independence p, i.e. destroying the order a 
bit faster (see Fig. [3J bottom row of Fig. 0] and most of the plots in Fig. [S] for further 
reference). 


4. Conclusions 


From physical point of view it is always an interesting question how details at the 
microscopic scale manifest at the macroscopic level. In the field of opinion dynamics 
such a macroscopic quantity is the opinion, defined in a case of binary models as the 
magnetization to. In this paper we examine six models that differ only in the way of 
selecting a group of influence but the size of this group remains fixed. Therefore there 
are no differences between models on the complete graph. For other topologies meth¬ 
ods for which range of interaction is shorter tend to show the lower critical value of 
the independence factor p below which to > 0 and above to = 0. Only two methods, 
RandBlock and Rand, give exactly the same results on all topologies and overlap MFA 
result found in [20]. Remaining four methods give different results and differences be¬ 
tween methods increase with the relative path length, i.e. are the most visible on the 
regular lattices. With decreasing relative path length the differences between methods 
vanish. One should notice that the average path length itself does not determine the 
differences between models, because results are virtually independent on the system size. 
What determines network properties is the relative path length defined as the ratio be¬ 
tween the average path length of considered network and the average path length of a 
random graph of the same size L and average degree (k). It should be noted that most 
of the real-world networks are characterized by relatively short paths 22] and therefore 
differences between models should be negligible. 
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We believe that our results contribute also to the discussion on differences between 
inflow and ouflow dynamics. As noted by Dietrich Stauffer: The crucial difference of the 
Sznajd model compared with voter or Ising models is that information flows outward: A 
site does not follow what the neighbors tell the site, but instead the site tries to convince 
the neighbors [4^- Of late, a debate on whether inflow dynamics is different from outflow 


dynamics has emerged 13, oei. Our findings indicate that not the direction of the 
information flow itself but the range of interactions is important, what coincides with 
results obtained by Castellano and Pastor-Satorras [17]. It is worth to notice that some 
of rules, investigated here, may be viewed as inflow and other as outflow dynamics. In 
particular, the NN method corresponds to the inflow dynamics. On the other hand, 
the Block method was inspired by the two dimensional and the Line model by the 
one dimensional outflow dynamics. Therefore both can be viewed as outflow dynamics. 
Both outflow rules ( Block and Line ) give the same results on scale-free and real Twitter 
networks, whereas the inflow rule NN gives lower value of the critical value of p. However, 
it seems that the critical value of p increases with the relative range of interactions and 
therefore it is understandable that NN (inflow) rule gives the lowest value of p. So 
perhaps one should not think about the direction (in or out) of the information flow itself 
but the range of interactions, what would coincide with results obtained by Castellano 
and Pastor-Satorras 113 . Summarizing, indeed inflow and outflow dynamics give different 
results but the reason is simply the difference in the range of interactions. 

As already mentioned in the introduction, the main motivation for this paper was the 
remark by Macy and Wilier that there was a little effort to provide analysis of how results 
differ depending on the model designs [hi . In the context of the problem posed here, 
it would seem that the structure of the group of influence may be important from the 
social point of view. However, as we have shown in the case of many complex networks, 
including BA and real networks, the importance of the group structure is often negligible. 

In this paper we considered only static networks (not changing in time), which is a 
common approach while studying dynamical processes like opinion spreading or diffusion 
of innovation. However, the characteristics of many real networks evolve in time and there 


is more and more data available on temporal networks |43j. If the changes take place 


at time scales comparable to those of studied processes, the temporal heterogeneities in 
such networks may lead to big differences in the dynamics of the processes, even if the 
networks appear similar from the static perspective [44^ . Thus, it could be worth to 
check what is the impact of the group structure in models put on top of real temporal 
networks. This issue will be addressed in one of the forthcoming papers. 
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